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ABSTRACT 

Around 2% of all A stars have photospheres depleted in refractory elements. This is hypothesized to arise from a preferential accretion 
of gas rather than dust, but the specific processes and the origin of the material - circum- or interstellar - are not known. The same 
depletion is seen in 30% of young, disk-hosting Herbig Ae/Be stars. We investigate whether the chemical peculiarity originates in a 
circumstellar disk. Using a sample of systems for which both the stellar abundances and the protoplanetary disk structure are known, 
we find that stars hosting warm, flaring group I disks typically have Fe, Mg and Si depletions of 0.5 dex compared to the solar-like 
abundances of stars hosting cold, flat group II disks. The volatile, C and O, abundances in both sets are identical. Group I disks are 
generally transitional, having radial cavities depleted in millimetre-sized dust grains, while those of group II are usually not. Thus we 
propose that the depletion of heavy elements emerges as Jupiter-like planets block the accretion of part of the dust, while gas continues 
to flow towards the central star. We calculate gas to dust ratios for the accreted material and find values consistent with models of disk 
clearing by planets. Our results suggest that giant planets of ~0.1 to 10M Jup are hiding in at least 30% of Herbig Ae/Be disks. 

Key words, accretion, accretion disks; protoplanetary disks; planet-disk interactions; stars: variables: T Tauri, Herbig Ae/Be; stars: 
abundances; stars: chemically peculiar 


1. Introduction 

Chemically peculiar A Boo stars have a solar-like photospheric 
abundance of volatile elements (C, N, O) but are underabundant 
by up to a factor of 100 in refractory elements (e.g., Fe, Mg, Si). 
The peculiarity affects up to 2 % of all late B through early F 
stars (e.g.. Gray & Corbally|[l998[ Paunzen||2001 >. Among the 
young Herbig Ae/Be stars, which host protoplanetary disks, the 
fraction is at least 33 % ( |Folsom et al.|201Z| . 

Stars of spectral type late B through early F have radiative 
envelopes, and consequently their photosphere is mixed with the 
interior on timescales ~ 1 Myr. This enables accreting material 
to dominate the photospheric composition for accretion rates as 


low as 10 13 M 0 yr 1 (Turcotte & Charbonneau]Jl993 


Turcotte 


2002) >. As the A Boo depletion pattern matches that of diffuse in¬ 
terstellar gas, it has been linked to the accretion of material from 
the circum- or interstellar medium, with radiation pressure pro¬ 
viding the mechanism that keeps dust from accreting while gas 
streams onto the star (Venn & Lambert 1990) Kamp & Paun- 
zen||2002 l Martinez-Galarza et al. 2009 . Since radiation pres¬ 


sure should affect all systems, it has been unclear why only a 
fraction of intermediate-mass stars display the A Bob peculiarity. 

Herbig Ae/Be stars are young B, A and F stars which host 
a protoplanetary disk and are accreting (Herbig 1960J | Waters | 
& Waelkens 1998). Based on their spectral energy distribu¬ 
tions, the disks are categorized as “group II” (cold, flat) and 
“group I” (warm, flaring; Meeus et a H2001| >. A further distinc¬ 
tion is made between disks showing emission from warm, small 
silicate grains (II, la) and those not showing such features (lb). 
Initially, group I was thought to evolve to group II as dust parti¬ 
cles grew and settled (Meeus et al.|200l) Dullemond & Dominik 


2004a b). The emerging view from recent debates is that group I 


disks are transitional, i.e., have radial cavities or gaps strongly 


depleted in millimetre-sized dust and perhaps gas (e.g., Keane| 


et al.|2014 jMaaskant et al. 2014 Banzatti & Pontoppidan|201 5[ 


Men u et al.|2015 van der Plas et al.|2015) l. As yet, no group II 

source is known to be a transition disk (TD), while all resolved 
TDs are group I sources. 

Cavities and gaps in transitional disks may arise due to dust 


being trapped in the outer disk by, e.g., giant planets (Whip 


pie 1972; Barge & Sommeria||l995; Rice et al. 

2006; 

Pinilla 

et al. 2012; Zhu et al. 2012), by gas and dust c 

earing due to 


photoevaporation ([Alexander et al. 2006a b] |201 


2011 20 


; Owen et aLj 


2|, or by the accumulation of gas and dust at the outer 


edge of a region with low ionization (Regalyet al.|201 2] Flock 


|et al.||20l5] >. The dominant mechanism should determine the 
volatile (gas) to refractory (dust) elemental ratios in the accretion 
columns that channel material onto the stellar photosphere. Most 
Herbig Ae/Be stars accrete at rates comparable to, or higher 


than, T Tauri stars (10 —10 ° M 0 yr , Garcia Lopez et al. 


2006 


Donehew & Brittain |201 1| Pogodin et al. |2012) >. We investi- 


gate whether stellar abundances correlate with disk properties, 
as might be expected from the above. 

2. Data 


Our sample, presented in Table A. 1 consists of the overlap be¬ 


tween the stellar abundance surveys and the compilation of disks 
described below. Metallicities or disk classifications for a few in¬ 
dividual sources were added from the literature, but in this regard 
the sample is certainly not complete. 

2.1. The elemental abundances 


We adopt the photospheric abundances from Acke & Waelkens 
( )20Q4[ > and Folsom et al. ( |2012j >. For overlapping sources, the 
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surveys are largely consistent, but we adopt the Folsom et al. 
values as the study is more consistent in the stellar parameter 
determination. For some Acke & Waelkens sources, we used val¬ 
ues based on the ionization state with the largest number of ana¬ 
lyzed lines. Some stars in the abundance studies host debris and 
not protoplanetary disks and were excluded. The abundances for 
FID 100546 are from Kama et al. (submitted) and for FID 34282 
from [Men'n et ak ( 2004 1. HD 101412 was identified as a A Boo 
star by Cowley et al. ( 2010] ), we use newer abundances from |Fol-] 
som et al. ( 201 2). Th e solar abundances shown in Fig.[T|are from 
Asplund et al. (|2009 1 . Values for the warm diffuse interstellar gas 
are from Savage & Sembach (19961. 


2.2. Properties of the disks 


Most of the disk group classifications are from Maaskant et al. 
( |2014j l. The references for radial disk cavities and gaps and for 
stellar accretion rates, spectral types and ages are given in Ta¬ 
ble [AT] HD 101412 is a TBoo star with a flat, cold dust disk 
consistent with a group II disk, but has a flaring, warm gas disk 
more reminiscent of group I systems (Fedel e et al.|2008 van der 
Plas et al.|2008| l. We have classified it as group I or transitional, 
as it is clearly not a pure group II source. HD 245185 is the only 
remaining massive dust disk in the 5 Myr old A Ori cluster and 
based on CO observations has a very low gas to dust mass ra¬ 
tio ( |Ansdell et al.|20f5) . HD 179218 is a distant (d = 244 pc) 
group I source which shows evidence for a characteristic dust 
emission radius larger than that of group II sources, although a 
cavity has not definitively been resolved (Leinert et al .|2004) l. 



Fig. 1 . Elemental abundances in the photospheres of Herbig stars host¬ 
ing group I (transitional/TD, open symbols) and group II (filled sym¬ 
bols) disks. The sun and the warm diffuse interstellar medium are shown 
in red. The full Herbig samples (circles) are shown, as well as their me¬ 
dians with sample standard deviations (diamonds with bars). 


3. Results 

Within our sample, all but two group I systems are confirmed to 
be transitional, with a large-scale reduction of the surface density 
of millimetre-size dust in the inner disk. As a working hypothe¬ 
sis, we consider the group I and transitional classifications to be 
identical (I/TD). We do not employ a quantitative cutoff between 
stars without and with A Boo abundance patterns, although some 
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of the stars in the sample display prominent A Boo anomalies. 
Rather, we separate the sources into group I and II disk hosts, 
and compare the stellar abundances between these subsets. 

Figure |T] shows the photospheric abundances of two volatile 
elements (C and O) and three rock-forming elements (Fe, Mg 
and Si) for the median and full sample of group I/TD and 
group II disk hosting stars. The solar photospheric abundances 
and those of the warm diffuse interstellar medium are given for 
comparison. The median elemental abundances are summarized 
in Table [T] The medians for the group II disk host stars are so¬ 
lar, while the medians for the group I/TD are solar-like for the 
volatiles, but depleted from solar by 0.5 dex for all refractory 
elements. The medians are unchanged if we exclude the Acke 
& Waelkens sample - all log 10 (X/H) values change by < 0.03, 
except Mg/H which increases by 0.13 for group II. 

Only two group II sources, HD 142666 and HD 144432, have 
strongly sub-solar refractory element abundances. Interferomet¬ 
ric observations reveal radial gaps on (sub-)au scales in both 


( |Chen et al.|2012||Schegerer et al. 201 3[ Menu et al.|2015| l. They 
are thus small-scale transitional disks with outer disks similar to 
group II sources. The group II/I unclear source HD 101412 - 
for which we adopted group I - has a narrow ring of CO emis¬ 
sion near the dust sublimation zone, also indicative of small- 
scale radial structure (Cowley et alTj|2012ji. The disks around 
HD 101412, HD 142666 and HD 144432 all have large radial 
variations of surface density at <1 au, with host stars depleted 
in refractory elements. As noted by Currie ( |2010) l, such systems 
may probe the earliest stages of disk clearing. 


Table 1 . Median log 10 (X/H) abundances and sample standard devia¬ 
tions (scatter within the subset) for group I and II disk host stars. 



C/H* 

O/H* 

Fe/H* 

Mg/H* 

Si/H* 

II 

-3.65 

-3.19 

-4.39 

-4.34 

-4.40 

cr 

0.07 

0.05 

0.17 

0.23 

0.24 

I/TD 

-3.57 

-3.14 

-5.04 

-4.84 

-4.95 

cr 

0.10 

0.10 

0.35 

0.35 

0.30 


3.1. The dust to gas ratio in the inner disk 

Group II disks are have no major depletion of dust or gas in their 
inner regions, while group I sources have surface density reduc¬ 
tion factors of 10 1 to > 10 6 for dust, and none to 10 for gas 
(|Bruderer et al.|20l4| JZhang et al.|2014| |Perez et al.|2015||van| 

der Marel et al.|2015) . A gap-opening giant planet reduces the 

surface density of large (millimetre-scale) dust by a larger factor 
than that of the gas ( |Rice et al.|2006||Zhu et al.|2012| l. This offers 
a natural explanation for a decreased photospheric abundance of 
refractory elements, as a large fraction of the dust mass can re¬ 
main trapped further out in the disk while gas and smaller grains 
continue to accrete. We now consider the photospheric compo¬ 
sition of Herbig stars as a proxy for the gas to dust ratio in the 
accreting material. The solar composition corresponds to a gas 
to dust ratio of 100. Noting the abundance of volatiles (C, O) as 
V/H and those of refractories (Fe, Mg, Si) as R/H, we express 
the A g /(j ratio in an inner disk feeding an accretion flow as 


A g/d = 100 x 


R/Ho 

R/H* 



(1) 


where V/H* and R/H* are the abundances in the host star. The 
results are shown in Fig.[2]for systems where at least one volatile 
and one refractory element has been measured. 
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related to the presence of an inner disk extending to the accre¬ 
tion radius (“pre-transitional” disks, see Espaillat et al. 2014j> 


6 8 10 
Age (Myr) 


12 14 16 


4. Discussion 

4.1. Evidence for giant planets 

The depletion of refractory elements in group I host stars is con¬ 
sistent with giant planets interacting with their disks. A planet 
with a mass <0.1 Mj up does not substantially alter the surface 
density of dust or gas, with A g /d » 100 at the planet position 
(e.g., iDong et al.|l2015). A planet in the mass range lM Jup < 
M p < 10 M j U p decreases the surface density of millimetre-sized 
dust more than that of the gas, with A g / d 


1000-10000 (Pinilla 


et al. 2012 Fig. 6). A more massive planet, M p > 10 - 15 M 


Jup» 

substantially decreases the surface densities of both the dust and 
the gas. Since the observations require a larger decrease of the 
dust than the gas in the accreted material, we propose that the 
A Bob-like Herbig stars host planets in the >0.1 to 10Mj up range, 
with the exact limits depending on the disk temperature and vis¬ 
cosity. Photoevaporation would deplete gas at least as efficiently 
as dust and is thus not consistent with the A Boo sources. We il¬ 
lustrate the typical group II and I inner disk structures and their 
relation to stellar photospheric abundances in Fig. [3] 

Fig.|T|shows that while the underabundance of refractory el¬ 
ements in group I/TD host stars is significant, not all are strong 
A Bob stars, and some have solar-like abundances. This may be 


and still containing enough mass in gas and dust to not be af¬ 
fected by the depletion of large grains in the gap containing 
the planet. Alternatively, such disks may host companions with 
masses > 10Mj up , which would block the accretion of both gas 
and dust from the outer disk. The stars in our sample have high 
accretion rates, which favours the massive inner disk scenario. 


Photosphere O 



Group II, full disk 

gas and dust flow onto star 



Fig. 2. The gas to dust ratios in the accreting material of group I (open 
circles) and group II (filled) disks, calculated from the stellar abun¬ 
dances with Eq.jT] 


3.1.1. HD 100546 

The highest inner disk gas to dust ratio in our sample, 1000, 
is found for HD 100546. We compare this photosphere-based 
value to the two-planet hydrodynamic disk model from Walsh 
et al. (2014 1 and Pinilla et al. ( |2015) >, and the physical-chemical 
disk model from Kama et al. (submitted). At 10 Myr, the age 
of the system, the hydrodynamic model yields A g;d ~ 1000 in¬ 
side 20 au. While this model was not configured to model the 
accretion zone, the gas to dust ratio of the inner cavity provides 
a zeroth-order approximation for the accreting material. The 
model from Kama et al. reproduces the spectral energy distribu¬ 
tion and a range of emission lines, and allows for A g /d = 1000 in 
the inner disk. The consistency of all three approaches suggests 
that the photospheric abundances reflect the inner disk material. 


Sun-like 


ISM-like 

(A Boo) 


Group I, trapping by planet 

gas (C, O) accretes 

dust (Fe, Mg, Si) trapped in disk 

Fig. 3. The stellar photospheric abundances and the inner disk gas and 
dust distribution for group I (transitional) and II sources. 


4.2. The occurrence rates of transitional disks, A Bod 
anomalies, and planets 

Based on a large homogeneous study of T Tauri stars with the 
SMA interferometer, | Andrews et ak[ ( |2011| ) found an occurrence 
rate of transitional disks of > 26%. As a first guess, we adopt 
this value for the Herbig Ae/Be systems. If the A Bob anomalies 
are due to the filtration of dust in transitional disks, we expect 
a similar occurrence rate of A Boo anomalies. Among the Her¬ 
big sample of Folsom et al. (2012), the TBoo occurrence rate is 
33 %, consistent with the estimated transitionality fraction. 

If the correlation of refractory depletion and disk structure 
is indeed due to giant planets, and if the transitional phase is 
short compared to the disk lifetime, the ~ 30% occurrence 
rate of transitionality and A Bob anomalies is a lower limit on 
the fraction of Herbig stars hosting giant planets in the 0.1- 
10 Mj up range. Current exoplanet surveys are not complete down 
to O.IMjup at Jupiter-like orbits. In the 3-14 Mj up mass range, 
the occurrence rate of planets at 5-320 au around A stars is 3 to 
30% at 95 % confidence (Vigan et al.|2012|, while the >4Mj up 
planet occurrence rate at 59-460 au tor B and A stars is <20 % 
(Niel sen et al.|2013] l. The radial velocity survey of Reffert et al. 
( 2015]> found an occurrence rate of ~30 % for planets >lMj up 
around intermediate-mass stars. This peaks towards super-solar 
metallicity ([Fe/H] > 0.2) and declines for stars more massive 
than ~2.5 M 0 . Allowing for undetected lower-mass planets, the 
occurrence rate of transitionality and A Boo anomalies is consis¬ 
tent with the giant planet population around A stars. The fact that 
the radial velocity planet occurrence rate peaks at a bulk stellar 
metallicity [Fe/H] > 0.2 suggests that A Boo-like depletions typ¬ 
ically decrease photospheric metallicities by more than 0.5 dex. 


4.3. The longevity of B, A, F star disks 

A large fraction of the disks studied in this paper are older than 
the canonical lifetime of 2-6 Myr for massive gas-dust disks 
(Haisch et al. 2001| >. Several sources are in the 10 Myr regime. 
At such ages, many systems already show low emission levels 
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consistent with debris disks (e.g., Hard y et al.|201 5| | Wyatt et al.| 
2015[>. This is particularly surprising since A stars are thought to 


lose their disks faster ( Ribas et al.|2015 1 and suggests that the gi¬ 
ant planet systems responsible for clearing large fractions of the 
inner disk may contribute to disk longevity by trapping material 
in the outer disk. It has been suggested that ~30 % of field stars 
retain their gas disks for for 10 Myr (Pfalzner et al. 2014j ). 

The frequency of the A Boo phenomenon among main se¬ 
quence B, A and F type stars is ~2% ( Gray & Corbally|1998[ >. 
These stars likely do not host massive protoplanetary disks and 
their A Bod peculiarity should disappear on a ~1 Myr timescale. 
However, a correlation has been found between a star show¬ 
ing A Boo peculiarities and hosting a debris disk ( jKing & Pat- 
ten |1992| Draper et al., submitted). This suggests that debris 
disks around early-type stars generally retain or produce large 
amounts of volatile-rich gas for up to Gyr timescales. We note 
that HR 8799, a 30 Myr old intermediate-mass star which hosts 
four giant planets of 5-7 Mj up at large orbits, as well as a multi- 
component debris disk, is a A Boo peculiar star (Marois et al. 


2008; |Su et al.||2009| |Marois et al.pOTOl |Murphy e t al.||2015| >. 

The photospheric abundances of early-type stars are a new win¬ 
dow into the long-term evolution of protoplanetary gas disks. 


5. Conclusions 

We present a correlation between the presence of a large radial 
cavity or gap in a protoplanetary disk and the depletion of refrac¬ 
tory elements in the photosphere of the host star. Such a deple¬ 
tion in intermediate-mass stars, known as the A Bob peculiarity, 
has been known for decades but has until now been unexplained. 
We propose that the trapping of large dust grains in the outer disk 
by giant planets provides a natural mechanism for increasing the 
gas to dust ratio of the accreted material, suggesting that as many 
as ~30 % of Herbig Ae/Be stars harbor giant planets. 

1. A typical late-B through early-F star hosting a group I or 
transitional disk is depleted in rock-forming elements (Fe, 
Mg, Si) by 0.5 dex compared to solar and group II host abun¬ 
dances. 

2. All group II disk hosts in our sample have solar-like stellar 
abundances. The two exceptional systems, HD 142666 and 
HD 144432, have small-scale radial dust gaps. 

3. The depletion of rock-forming elements in stars hosting 
group I disks is likely due to giant planets which trap large 
dust grains in the outer disk. The diversity of disk structures 
and planet masses explains why A Boo-like stars have a range 
of depletions, rather than clustering at particular values. 

4. A deficiency of rock-forming elements in the photosphere of 
a disk-hosting B, A or F star implies the presence of a radial 
dust cavity or gap, and a giant planet. 

5. Most as-yet uncharacterized central stars of transitional Her¬ 
big Ae/Be disk systems will likely be found to have a low 
abundance of rock-forming elements. 
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Table A.l. The systems studied in this paper. 
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